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Line-of-sight path delay calibration accuracies of 1 mm are needed to improve 
both angular and Doppler tracking capabilities. Fluctuations in the re fr activity of 
tropospheric water vapor limit the present accuracies to about 1 nrad for the angu- 
lar position and to a delay rate of 3 x 10" 13 sec/sec over a 100-sec time interval for 
Doppler tracking. This article describes progress in evaluating the limitations of the 
technique of water vapor radiometry at the 1-mm level. The two effects evaluated 
here are: (1) errors arising from tip-curve calibration of WVR’s in the presence 
of tropospheric fluctuations and (2) errors due to the use of nonzero beamwidths 
for water vapor radiometer (WVR) horns. The error caused by tropospheric wa- 
ter vapor fluctuations during instrument calibration from a single tip curve is 0.26 
percent in the estimated gain for a tip-curve duration of several minutes or less. 
This gain error causes a 3-mm bias and a 1-mm scale factor error in the estimated 
path delay at a 10-deg elevation per 1 g/cm 2 of zenith water vapor column density 
present in the troposphere during the astrometric observation. The error caused by 
WVR beam averaging of tropospheric fluctuations is 3 mm at a 10-deg elevation per 
1 g/cm 2 of zenith water vapor (and is proportionally higher for higher water vapor 
content) for current TO beamwidths (full width at half maximum of approxi- 
mately 6 deg). This is a stochastic error (which cannot be calibrated) and which 
can be reduced to about half of its instantaneous value by time averaging the radio 
signal over several minutes. The results presented here suggest two improvements 
to WVR design: First, the gain of the instruments should be stabilized to 4 parts 
in l(ft over a calibration period lasting 5 hours, and second, the WVR antenna 
bea mwidth should be reduced to about 0.2 deg. This will reduce the error induced 
by water vapor fluctuations in the estimated path delays to less than 1 mm for the 
elevation range from zenith to 6 deg for most observation weather conditions. 
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I- Introduction 

Future missions will benefit from spacecraft tracking 
with 100-prad accuracy in an angular position. The 
Cassini radio science team requires a delay-rate accuracy 
of 5 x 10“ 16 sec/sec over 1000 sec for gravitational wave 
searches using Doppler tracking. Fluctuations in wet tro- 
pospheric refractivity limit the present tracking capabili- 
ties to about 1 nrad for angular position and 5 x 10“ 14 
sec/sec over 1000 sec for Doppler tracking, which corre- 
sponds to a 0.1-mm/sec uncertainty in spacecraft radial 
velocity. 

Angular tracking is done with very long baseline inter- 
ferometry (VLBI), a technique that measures the differen- 
tial phase between two DSN antennas of an electromag- 
netic signal originating from a radio source. By relating 
this measured phase difference to geometrical path delays, 
astrometric parameters can be estimated. Inhomogeneities 
in tropospheric refractivity cause unmodeled errors in the 
path delay at about the 1-cm level for a path delay at 
zenith over a period of several hours. This error t6r- 
rupts angular position estimates at about the 1-nrad level. 
To achieve 100-prad angular position accuracy, fluctua- 
tions in path delays must be calibrated at the 1-mm level 
[ 1 , 2 ]. 

Charged particles (both in the Earth’s ionosphere and 
in solar wind) are the main source of error for spacecraft 
gravitational wave searches using Doppler tracking at S- 
band (2.3 GHz) (all missions prior to Galileo). Searches 
planned for Galileo at X-band (8.4 GHz) will be limited by 
fluctuations in the plasma and the troposphere at about 
the same level, a delay rate of approximately 5 x 10“ 15 
sec /sec over 1000 sec [3], In searches planned for Cassini 
at K-band (32 GHz), the increased observational frequency 
will reduce the plasma-induced error to approximately 
5 x 10“ 16 sec/sec over 1000 sec. To take advantage of 
this increased sensitivity, tropospheric fluctuations must 
be calibrated at the submillimeter level. 

Water vapor radiometers (WVR’s) have been suggested 
for measuring line-of-sight path delays due to tropospheric 
water vapor. WVR’s work on the principle that the radio 
power collected by a WVR antenna is proportional to the 
brightness temperature of the sky in the antenna point- 
ing direction. Path-delay retrieval algorithms relate the 
brightness temperatures measured at two (or more) fre- 
quencies near the 22.6-GIIz water vapor absorption line 
to a path delay associated with tropospheric water vapor 
along a line of sight in the same direction [4]. Brightness 
temperatures at two or more frequencies are required to 


subtract the effect of liquid water absorption from the to- 
tal absorption. Unlike water vapor, liquid water does not 
affect the refractivity at microwave frequencies. 

The ability of WVR’s to determine the absolute path 
delays, or to track path delay changes induced by tropo- 
spheric fluctuations, has been tested recently with mixed 
results. 1 3 The accuracy of WVR calibration using radio- 
sondes is limited at a level of approximately 10 percent 
(which causes a 0.6-cm error for a 6-cm zenith path delay) 
due to uncertainties in both the radiosonde data and path 
delay retrieval algorithms [4]. 2 Studies that calibrated 
VLBI time series with differenced WVR delays using co- 
pointed antennas reported reduced rms residual delays at 
high elevations, 3 but the residual delays actually increased 
(up to 20 psec, corresponding to a 0,6-cm path delay er- 
ror) at elevations below 50 deg. These results strongly 
suggest that to meet the future mission requirements, 
the character of various error sources that affect the ac- 
curacy of water vapor radiometry must be much better 
understood. 

The accuracy with which WVR’s estimate tropospheric 
path delays is determined by error sources that include 
(but are not restricted to) inaccuracies in calibration of the 
WVR gain, uncertainties in the path-delay retrieval algo 
rithms and atmospheric absorption modeling, radiometer 
noise, the effect of the WVR location relative to the ra- 
dio telescope, and the beam intensity distribution. The 
error sources will be discussed later. The errors can origi- 
nate in the instrument or in the atmosphere and can cause 
bias, scale, or random errors in the retrieved path delays. 
The effect of a constant bias on path delays for angular 
tracking can be eliminated by differencing between obser- 
vations. Biases have no effect on delay rates used for grav- 
itational wave searches. Unmodeled variability in spatial 
distributions of atmospheric parameters together with the 
system noise have been recognized as the principal noise 
mechanisms that limit the ability of WVR’s to monitor 
tropospheric fluctuations. Only recently has it been recog- 
nized that for any realistic WVR design, the tropospheric 
dynamics will also influence the WVR accuracy. 


T. J. Vesperini, M Staplelon WVR Experiment-Part I: Results," 
JPL Interoffice Memorandum (internal document), Jet Propulsion 
Laboratory, Pasadena, California, March T, 1988. 

^ S. Keihm, “Water Vapor Radiometer Intercomparison Experiment: 
Platteville, Colorado, March 1-14, 1991," Final report, JPL Task 
Plan 30-3289 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, July 1991. 

3 C. Edwards, “Water Vapor Radiometer Line-of-Sight Calibration 
Capabilities," JPL Interoffice Memorandum 335.1-90-015 (inter- 
nal document), Jet Propulsion Laboratory, Pasadena, California, 
March 30, 1990. 
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The main goal of this article is to investigate how the 
tropospheric dynamics affect the WVR’s ability to track 
tropospheric fluctuations. Because WVR’s are imperfect 
instruments, one of the goals of this article is to quantify 
relations between WVR design parameters and path-delay 
retrieval accuracy. To understand how the errors affect 
path delay estimates, an analysis of the effect of WVR 
error sources on path delay retrieval is presented. The 
article then focuses on two specific errors that arise from 
tropospheric fluctuations: errors in WVR instrument gain 
calibration from tip curves and errors caused by WVR an- 
tenna beam averaging. Tropospheric fluctuations cause 
errors in the estimated gain because data at different tip 
directions are analyzed using mapping relationships valid 
for a temporally constant and spatially homogeneous tro- 
posphere. Tropospheric fluctuations cause unmodeled de- 
partures from this picture, which induces errors in the es- 
timated gain and ultimately in the estimated path delays. 
WVR antenna beam averaging causes errors in the esti- 
mated brightness temperature in the direction of a radio 
telescope because data recorded by the WVR’s are beam 
averaged around the WVR pointing direction (full beam 
widths of the present state-of-the-art WVR’s are between 
4 and 10 deg), whereas radio telescopes measure the tro- 
pospheric effects along the line of sight to a distant radio 
source. The averaging corrupts the accuracy of the esti- 
mated brightness temperature in the radio source direction 
for all realistic WVR designs. 

The article is organized as follows: Section II dis- 
cusses how various error sources affect path delay esti- 
mates. Some basic equations describing the conversion of 
the recorded data to the line-of-sight path delay are given. 
Section III calculates the error in the WVR measurement 
incurred by using tip curves in an inhomogeneous tropo- 
sphere. This error will exist as long as the WVR’s are cal- 
ibrated by using tip curves. Section IV studies the effect 
of beam averaging on WVR measurements for collocated 
antennas. The effects of antenna copointing and beam in- 
tensity distribution are discussed. The aim here was not 
to derive exact numbers for any specific WVR design but 
rather to provide rough error estimates. Section V is a 
summary with recommendations for WVR gain stability 
requirements and antenna beam width design to comply 
with the 1-mm path-delay accuracy requirement, as well 
as plans for further studies. 

II. WVR Error Sources and Their Effect 
on Path Delay Retrieval 

To put the above-mentioned effects in perspective, this 
section gives an overview of how various error sources af- 
fect the estimated path delay. The conclusions of this 


overview are summarized in Fig. 1. Error sources result 
from inaccuracies in WVR measurements or data inter- 
pretation and can cause bias, scale, or random effects in 
the retrieved path delay. The effect of a constant bias or 
linear trend on the tracking of changes in the path delay 
can be eliminated with astrometric parameter estimation. 
A scale error is an error that is directly proportional to 
the path delay. It is unclear at the present time to what 
extent radiometric data reduction can filter out the effect 
of a systematic scale error on astrometric parameter esti- 
mates. Stochastic errors cannot be reduced by parameter 
estimation. The effect of error sources on the estimated 
path delays can be analyzed by using a path delay retrieval 
algorithm that relates the line-of-sight path delay ( L v ) due 
to water vapor in the troposphere to the brightness tem- 
perature ( Tb ) in the same direction. For example, Eq. 
(20) of [4], 

L v — ao + QiTb,i + a 2 T B , 2 (1) 

expresses L v asa linear combination of the brightness tem- 
peratures Tbj (j = 1,2) at frequencies /i = 20.6 GHz 
and / 2 = 31.4 GHz. In Eq. (1), a/s are coefficients that 
have been determined by a regression analysis of the WVR 
data [4]. The value of ao quantifies the effect of the dry 
component of the atmospheric absorption; a\ quantifies 
the effect of water vapor absorption; and the term a 2 T B>2 
(a 2 ^ —cii (fi/fi) 2 ) subtracts the contribution due to tro- 
pospheric liquid water from the total absorption (liquid 
water causes a negligible path delay). For L v in centime- 
ters and T bj in kelvins, typical regression coefficients are 
ao ~ —1.6, ai ~ 0.66, and a 2 a —0.3. Typical Tbj ’ s 
are between 10 and 100 K. The value of L v scales as 
it, (cm) ~ 6 N v f sin E ) where N v is the water vapor col- 
umn density (in g/cm 2 ) at zenith, and E is the observed 
elevation. For typical N v between 1 g/cm 2 and 4 g/cm 2 , 
the zenith path delay ( L V}Z ) is between 6 and 24 cm. 

Taking the differential of Eq. (1) determines the error 
in the retrieved path delay, SL V) in terms of errors in the 
regression coefficients, tfa^’s, and the errors in Tbj 7 s as 

2 

6L V = 6ao + ^2f{^ a jTsj + ajSTsj) (2) 
;=i 

where the subindex j = 1,2 refers to the WVR frequency 
channels. This article focuses primarily on errors in Tbj 
induced by tropospheric fluctuations. However, for the 
completeness of the discussion and because they are so 
large, sources of errors in a/s are briefly summarized first. 
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The biggest errors in aj 5 s come from two sources: (1) inac- 
curate modeling of the absorptivity of water vapor, a v [i.e., 
from errors of the dependence of ot v — a v (T,p, /?, /) on the 
ambient temperature T, pressure p, water vapor density p, 
and frequency /], and (2) uncertainties in the spatial and 
temporal distributions of atmospheric pressure, tempera- 
ture, and water vapor and liquid present in the troposphere 
during the observation along the line of sight. Errors in the 
aj's cause scale errors in L V} of the type SL V cz L v Sa\/ai. 
The error caused by inaccurate modeling of ot v (T, p, p, /) 
is systematic, about 10 percent for the current absorptiv- 
ity models. Atmospheric profile uncertainties depend on 
season, site, weather, time of day, and line of sfght. The 
uncertainties are difficult to model and are the reason why 
there is no one-to-one correspondence between the bright- 
ness temperature and the path delay. The error caused 
by uncertainties in atmospheric profiles is between 2 and 
4 percent, depending upon the specific retrieval algorithm 
used. Several possibilities pertaining to the feasibility of 
reducing the 8a j ' s will be noted in Section V. 


The error in the brightness temperature, 8T B j , comes 
from two sources: (1) the error STaj in the measured 
WVR antenna temperature Taj> and (2) an error in the 
interpretation of Taj. The latter contribution to ST B j 
originates in the fact that WVR’s do not measure T B j di- 
rectly, but rather they record a signal, Vaj (from which 
Taj is extracted), and the obtained Taj is used to esti- 
mate the brightness temperature T B j. The error 8Taj 
in the WVR-measured Taj can be expressed in terms of 
the WVR parameters, as follows: All WVR’s use an in- 
ternal reference, such as blackbody radiation or a noise 
diode, to enable a subtraction of the contribution of the 
system temperature from the recorded data, The recorded 
Vaj is proportional to the difference between Taj and the 
reference load temperature, T rt j, 

Vaj = g(T AJ -T rtJ ) (3) 


where g is the WVR gain, Inverting Eq. (3) to obtain Taj 
and taking a differential of the resulting equation yields the 
error STaj in the measured Taj 


8Va - 

&Taj = -f*- + 6T reJ - ^V A j 
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ence load temperature (6T re j) and the WVR gain (<$<7). 
Figure 1 shows that for the path delay error to be less 
than 1 mm, the noise (or more generally, any unmonitored 
drifts) in all temperature-like quantities must be less than 
0.2 K. In practice, the effect of the system noise on path 
delay estimates can be reduced (at the expense of time res- 
olution) by increasing the signal integration time. As long 
as T re f remains constant, the path delay error caused by 
using an incorrect value of T re j is a constant bias (whose 
eflect on astrometric estimates can be eliminated by dif- 
ferencing between the observations). 

Because the antenna voltage Vaj oc Taj — T re j t the 
error caused by Sg [the last term on the right-hand side 
of Eq. (4)J consists of two terms. The first is oc T rt j , the 
second is oc Taj . Note that since the typical T re j ~ 300 K 
is bigger than Taj (Taj ~ T B j 10 to 100 K), the pres- 
ence of Trtj in Vaj enhances the effect of the gain $f?5r 
(whatever its origin may be) on 6L V . By using Eq. (4) 
in Eq. (2), the first term (i.e., the term oc T re f) leads 
to 6L V oc a, T re j 6g/g. For stable gain (i.e., for a con- 
stant difference between the estimated and the WVR true 
gains), this is a constant bias. For unmonitored gain fluc- 
tuations, this is a random error, which, in order to satisfy 
the i-mm path delay requirement (see Fig. 1), must be 
6 g/g < 0.08 percent. The second term (i.e., the term 
oc Taj) is a scale error, 6L V ~ L v 6g/g , which is system- 
atic for a stable gain and time varying for an unstable gain. 
This error depends on the tropospheric humidity and the 
observed elevation. For the error to be less than 1 mm, the 
gain error must be Sg/g < 0.1 sin E/L v>z (where L v z is 
the zenith path delay in centimeters and T VtZ = 6 cm for a 
troposphere with 1 g/cm 2 of water vapor column density 
at zenith). Figure 1 shows that for the scale error to be 
less than 1 mm at a 10-deg elevation when = 6 cm, 
6g/g must be less than 0.26 percent. For the error to be 
less than 1 mm at E = 6 deg when L VfZ = 24 cm, 6g/g 
must be less than 0.04 percent. 

The gain error can originate from several sources. The 
most troublesome of these are unmonitored instrumental 
drifts on a time scale of 100 to 1000 sec 4,5 (stochastic 
fluctuations on time scales much shorter than the radio 
observations can be incorporated into the system noise). 
However, since they originate in the instrument, the drifts 
should be controllable by improved instrument stabiliza- 


= + 6T rt j + —(T re f - T A j) (4) 

9 9 

Thus, STaj comes from three WVR parameters, the 
system noise (modeled as the uncertainty in the equivalent 
temperature, SVAj/g), and the uncertainties in the refer- 


4 G. M, Resch (Tracking Systems and Applications Section) and 

S. Keihm (Microwave Observational Systems Section), personal 
communication, Jet Propulsion Laboratory, Pasadena, CaEfbrnia, 
19 81 . . ; ; 

5 G. Parks, C. Ruf, and S. Keihm, “Advanced Water Vapor Ra- 
diometer: Definition Phase Study” (internal document), Jet Pro- 
pulsion Laboratory, Pasadena, California, November 29, 1990. 
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tion in advanced WVR designs. Tropospheric fluctuations 
are another source of 6g for WVR’s calibrated using the 
tip curves. Even though this 6g is stochastic in origin, the 
ensuing difference between the estimated and the actual 
WVR gains causes a systematic (bias and scale) error in 
L v . The error depends on tropospheric humidity and cal- 
ibration strategy and will be calculated in Section III of 
this article. 

The other contribution to 6T B j is the error made in in- 
ferring Tbj from Taj- The simplest and most often used 
relationship between Tbj and Taj is that Tbj = Taj - 
This relationship can be in error because of inaccurate an- 
tenna pointing and spatial separation between the WVR 
antenna and the radio telescope. Another source of uncer- 
tainty is the effect of averaging tropospheric fluctuations 
over WVR beam intensity distribution. Section IV de- 
rives an expression for the error in the estimated Tbj that 
was caused by the beam averaging of tropospheric fluctu- 
ations for collocated and copointed antennas. The error 
will increase with decreasing elevation more rapidly than 
L V} i.e., faster than a simple scale error. Since, as a result 
of this increase, low-elevation data will be weighted more 
heavily than they should be during VLBI data reduction, 
astrometric parameter estimates will be impacted. 


III. Error In the WVR Gain Estimated From 
Tip Curves Due to Tropospheric 
Fluctuations 

This section presents a calculation of the error in the es- 
timated gain (£) of WVR’s induced by tropospheric fluctu- 
ations during the WVR calibration using the tip curves. 6,7 
Tip curves use the elevation dependence of the sky bright- 
ness temperature [. Tb t i = T B (Ei), where Ei is the tip 
elevation] to calibrate the WVR gain. The WVR gain 
is determined by fitting the WVR recorded signal Vaj 
(V A t i = V A (Ei)). If the elevation dependence of T B j> 
and therefore of Vaj, were known, this calibration proce- 
dure would be limited only by thermal measurement noise. 
Spatial and temporal fluctuations of atmospheric water va- 
por cause the actual distribution of water vapor to depart 
from the static distribution assumed in fitting the Va/ s. 


® J. Z. Wilcox, “The Standard Deviation of WVR Gain Estimated 
from Tip Curves due to Wet Troposphere Fluctuations,” JPL 
Interoffice Memorandum 335.6-91-032 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, December 19, 1991. 

7 J. Z. Wilcox, “The Difference Between Two Successive WVR Gain 
Estimates FVom Tip Curves due to Wet Troposphere Fluctuations,” 
JPL Interoffice Memorandum 335.6-91-033 (internal document), 
Jet Propulsion Laboratory, Pasadena, California, December 20, 
1991. 


Therefore, an error due to tropospheric fluctuations is in- 
troduced into the gains estimated from the tip curve data. 
All delays calculated with the derived gain will, therefore, 
also be in error. In this section, the covariance of the gain 
estimates is determined in terms of the covariance of the 
tropospheric opacity. This is done by performing a tip- 
curve analysis of modeled data and evaluating the opacity 
covariance by using the Kolmogorov turbulence model [1] 
for the wet troposphere. 

A. Tip-Curve Analysis 

The tip-curve data were modeled by using the optically 
thin tropospheric approximation for the standard radia- 
tion transport equation, neglecting the effect of the Earth’s 
curvature, ray bending, and nonzero WVR beamwidth, 
assuming negligible time elapsed during each tip curve 
sequence, and neglecting the time variation of all other 
model parameters except the tropospheric opacity. The 
recorded signal Va,*( 0 for tip curve epoch t at elevation 
Ei is then expressed as [5] 

V A ,■(<)= 9 {T B j{t)-T r . s ) 

= 9 (Tc e" T ‘« + T m ( 1 - e~ T 'W) - T rtJ ) 

— 9 (Tc + Tmc Ti(t) — T re /) (5) 

where the subscript i refers to ith elevation, g is the 
WVR gain, T reJ (T rtJ ~ 300 K) is the reference tem- 
perature discussed in Section II, Tc ~ 2.8 K is the cos- 
mic background temperature, Tmc — Tm — Tc where 
Tm ~ 270 to 280 K is the average atmospheric temper- 
ature [5], and r*(*) (tv(<) = r(Ei 1 t)) is a time-varying line 
of sight opacity at elevation E{. The linear approximation 
[the second expression on the right-hand side of Eq. (5)] 
is a good approximation for most optically thin (r* <0.5) 
tropospheres of interest, with the added benefit of mathe- 
matical simplicity. 

In the standard tip curve analysis, the method of least 
squares [6] is used to obtain the gain estimate ( g ) in terms 
of the tip data. Appendix A discusses the tip curve fitting 
in detail. The data are fit to a static (temporally averaged) 
version of Eq. (5) by using the mapping function (r,) = 
r z A{ ) where (...) designates the statistical average, r 2 is 
the averaged opacity mapped to zenith, and the air mass 
Ai = 1 / sin E { . For N elevations, the fitting leads to an 
equation of the following type; 

N 

m = E c * ( 6 ) 

i=i 
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where c,-*s are coefficients that depend on the assumed Tc, 
7 mc> and T re j f and on the so-called variance-covariance 
matrix W~ 1 = cov(V^(<), V^t')) (^ ee l®3 ^ Appendix 
A). Taking Eq. (6) at t and t f and substituting Eq. (5), the 
covariance of the gain estimates separated by the time in- 
terval T = t—t l is obtained in terms of opacity correlations 
as 

. — N 

co v(g(t),g{t')) = ^2 Ci Cj cov(Va,,(0Va, .(<')) 
ij = 1 

- 9 2 The Y2 °i c ) cov (ri(0^(<')) (7) 

•j=i 

The opacity covariance cov(r,(f)rj(< / )) was evaluated by 
neglecting fluctuations in the dry component of r»(f) (dry 
troposphere contributes less than about 30 percent of the 
total opacity fluctuations) 5 and by describing the fluctifk- 
tions in the wet component by Kolmogorov turbulence [1]. 
Specifically, the wet contribution to r,- was expressed as 
the line-of-sight integral = f* thv a v (rj,f) dr,*, where 
o v (r,t) is the tropospheric absorptivity due to water va- 
por per unit length at r, dr, = A{dz is the path increment 
along the line of sight at elevation Ei> and h v is the wet tro- 
posphere height. Using a v (r t t) ~ x(r>t) T v t z/L VfZ} where 
X — index of refraction — 1, and t V jZ and L VtZ are wet 
opacity and path delay at zenith, respectively, Appendix 
B shows explicitly the integral expressions that relate 
cov(ri(f)rj(f / )) to the structure functions for x(^0- (Note 
that L V}Z ~ 6 cm, and r VfX ~ 0.04 and 0.02 per 1 g/cm 2 
of zenith water vapor column density, and the dry opacity 
Td fZ — 0.017 and 0.04, at 20.6 GHz and 31.4 GHz, respec- 
tively.) The refractivity structure functions were evaluated 
by generalizing the Kolmogorov turbulence expression [1] 

for the structure function ^(x(r) — x(^+ R)) 2 ^ to inho- 
mogeneities correlated both spatially and temporally [l]: 9 

D X (R, T) = / ( X (r, t) - X (r+ R, t + T , )/) 

_ NZ C 2 I R+II^l (8) 

1 + (\M + vT\/L,)W \ } 

where R and T are the spatial and temporal intervals over 
which the structure function is evaluated and v is the wind 
velocity. T*he role of Eq. (8) in VLBI data reduction was 


8 G. E. Lanyi, personal communication, Tracking Systems and Ap- 

plications Section, Jet Propulsion Laboratory, Pasadena, Califor- 

nia, 1991. 

8 See Footnote 7. 


discussed in [1]. By using the standard deviation for re- 
trieved path delays at average DSN conditions, the turbu- 
lence strength was shown to be C = 2.4 x 10“ 7 m -1 ^ 3 for 
a tropospheric slab with a water vapor column density of 
N v ~ 1 g/cm 2 at zenith (corresponding to approximately 
6 cm of wet path delay) and height h v = 1 km [1], For 
h v = 2 km, the recalculated C = 1.1 x 10 _7 m _1 / 3 , which 
is the value used in this article. 10 The turbulence satu- 
ration scale length was taken to be L $ — 3,000 km [1]. 
The temporal correlation depends on the wind velocity v 
(Reference [1] shows that if one identifies T with \R\/v t 
the spatial correlation between two parallel lines of sight 
separated by the distance |/q projected on the Earth’s sur- 
face is equal to the temporal correlation of a single line of 
sight at time t and later t+T.) This article used v = 10 
m/sec, which is a typical wind speed at the Goldstone DSN 
antenna site. 

B. The Estimated Gain Error 

The standard deviation of the estimated gain, <7^, is 
equal to the square root of the covariance given by Eq. 
(7) for t = For t ^ t f y Eq. (7) gives the gain covariance 
as a function of the time interval T = \t — f] between two 
successive gain estimates, (r${T) = (cov (g(t) } g(t + T))) x ^ f 
Before discussing the numerical results, note that Eq. (7) 
depends on the number ( N ) of tip elevations. Evaluating 
Eq. (7) for N = 2, 3, and 4, the covariance was found to 
depend on the tip range and be nearly independent of the 
elevation distribution within the range. 11 That is, the co- 
variance is determined by the least-correlated tropospheric 
inhomogeneities (i.e., by correlations between the lines of 
sight associated with the minimum, E m i n) and maximum, 
E t naan tip elevations). In what follows, the errors will be 
shown versus E m i n (with E max at zenith) per 1 g/cm 2 of 
zenith column density of water vapor present in the tropo- 
sphere during WVR calibration. Note that because they 
are proportional to N Vy the errors are minimized by cali- 
bration in dry (and stable) weather. 

The value of is shown in Fig. 2 as a function of 
E m i n * The error has a flat minimum of approximately 
0.26 percent between approximately E m i n = 10 deg and 
30 deg. Below a 10-deg elevation, the error increases with 
decreasing E m i n because the decorrelation between the 
tropospheric fluctuations associated with E m i n and 90-deg 
lines of sight increases more rapidly than does the air mass 
difference. Above a 30-deg elevation, the opposite is true. 
It can be shown that as E m i n approaches 90 deg, the er- 
ror increases as (90 — i?mm)“ 3 ^ 2 (which is slower than it 


10 See Footnote 6. 

11 See Footnote 6. 
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would be if the fluctuations were completely random, in 
which case the error would increase as (90 — E m in) )• 
Thus, to minimize the gain error, Z7 m ,* n should be between 
10 and 30 deg. 

The least-squares fit coefficients c^s depend on the 
variance-covariance matrix W~ l , In the so-called con- 
sider analysis, the observable errors are assumed to be 
uncorrelated, and W~ l is approximated by a unit diag- 
onal [1]. Using the observable variance-covariance matrix 
W~ l minimizes the variance of the estimated gain, Eq. (7) 
(i.e., it minimizes <Tg). The errors calculated by using the 
unit and observable variance-covariance matrix W~ l are 
shown as solid and broken line curves, respectively, in Fig. 
2. Note that the two curves are practically the same in 
the minimum region and differ by 8 percent at most in 
the wings. This indicates that using the full covariance- 
variance matrix W~ l does not significantly improve the 
estimated gain accuracy. 

Figure 3 shows the temporal development of 0g{T). 
The value of <Tg(T) decreases at T « T corr approximately 
as <?g(T) ~ Gg\J 1 - (2T /Tcorr), where T corr is the decor- 
relation time. (Note that in Fig. 3, crg(T corr /2) ~ <?gf2.) 
The value of 'Tcorr depends on v and Emin as Tcorr — 
h v /(v tan E m in)' That is, T cor r is the time it takes for a 
“frozen” troposphere [1] to pass through the tip range be- 
tween E m in and zenith. For v = 10 m/sec, and E m in — 30 
deg, T C orr is about 7 min. The single-gain estimate errors 
(tf^’s) become independent from each other when the time 
T between subsequent tip curves exceeds T corr . Therefore, 
if one wishes to minimize the estimated gain error by tip 
curve repetition, the tip curves should be separated by a 
time interval greater than 7 min. 

C. The Estimated Path Delay Error 

It has already been discussed in Section II that a WVR 
gain error induces two types of errors in the estimated 
L v : a bias and a scale error. Using Eq. (4) in Eq. (2) 
and designating the bias and scale errors as A L v and SL V , 
respectively, the two errors are 

AL v (cm) ~ (ai + a 2 ) T re / 

~ 120— ~ 0.3 N v cat (9a) 

9 

6L v (cm) ~ (a, Ta,i + ct 2 Ta, 

~ L v — ~ 0.016 ( 9 b) 

g sm E 


where the gain estimates at 20 6 GHz and 31.4 GHz were 
assumed to be correlated (i.e., the gains in the two WVR 
channels were determined during the same tip sequence) 
and their magnitudes the same. T re j = 300 K, N VtCa i and 
N Vi t oi s are the number of grams per cm 2 of the column 
density of water vapor at zenith during the WVR calibra- 
tion and radio observation, respectively, E is the elevation 
of the observation, and the last expressions on the right- 
hand sides correspond to the minimum <Tg/g ~ 0.26 per- 
cent. The A L v corresponding to N v>ca i = 1 g/cm 2 (i.e., 
to L VtJ = 6 cm) is shown in Fig. 2. The minimum A L v 
is approximately 3 mm. Note that as long as the WVR 
gain remains constant, A L v is also constant, which makes 
it possible to remove its effect on astrometric estimates 
and delay rates by differencing between observations. For 
unstable gain, the gain changes must be monitored (to 
achieve a 1-mm path delay accuracy) with 0.08 percent ac- 
curacy. The dependence of the scale error (strictly speak- 
ing, 6L V will be a pure scale error only in optically thin 
tropospheres) on elevation was shown in Fig. 1. The value 
of 6L V caused by a 0. 26-percent gain error at a 10-deg el- 
evation when N Vi0 i„ = 1 g/cm 2 is approximately 1 mm; 
when N Vi ob s = 4 g/cm 2 , the error is 4 mm. To reduce 6L V 
to 1 mm, the gain error will have to be reduced by us- 
ing a different calibration technique, tip curve repetition, 
or a parameter estimation during the data analysis. How 
this can be accomplished for a stable WVR is discussed in 
Section V of this article. 

A note should be made here on the dependence of A L v 
and 6L V on T re j . From Eq. (9a), it would seem that A L v oc 
Tre/. However, by performing the least-squares analysis, 
one finds that the tropospheric fluctuation-induced <jg is 
oc 1 /T re j . This cancels the dependence of A L v on T re j and 
makes 6L V oc 1 /T rt f- [Note, however, that A L v caused by 
an instrumental gain drift will be a T re f , as given by the 
first expression on the right-hand side of Eq. (9a).] 


IV. Error In the Estimated Brightness Due 
to the WVR Nonzero Beamwidth^ 

Retrieval algorithms relate the line of sight L v to the 
brightness temperature ( Tq ) in the same direction, where- 
as data recorded by the WVR’s axe beam averaged around 
the WVR pointing direction. For collocated and copointed 
WVR antennas and radio telescopes (the telescope points 
along the line of sight to the radio source), the copointing 
introduces two types of errors into the estimated 7#: a 
systematic error due to the nonlinear dependence of air 
mass on elevation and a random error due to WVR beam 
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averaging of tropospheric fluctuations. 12 The systematic 
error can be calculated for known beam intensity distribu- 
tions (beam shapes) and water vapor content in the tro- 
posphere. Its effect on path delay estimates can be elimi- 
nated by pointing the WVR to a slightly higher elevation 
so that the radio source lies in the direction of the centroid 
of the distribution of WVR beam brightness. In this sec- 
tion, the direction of the brightness centroid is calculated 
by WVR beam averaging of the air mass [see Eq^ (12) 
for the centroid definition}, and the random error is deter- 
mined by using the tropospheric opacity statistically as in 
Section III. It has been suggested that to simplify WVR 
beam steering for collocated antennas, radio telescope and 
WVR antennas should be copointed. 13 After correcting 
the WVR data by subtracting from them the systematic 
error, the “corrected” data would be used to estimate path 
delays in the direction of the beam’s geometrical center. 
Therefore, the random error has been evaluated also for 
this geometrical center pointing case (and found that it 
is bigger than the random error for the brightness cen- 
troid pointing, by an amount that depends on the WVR 
beamwidth and elevation). 


tively, and B c (E y <p) is the WVR antenna beam radiation 
pattern 


B c (E,<p) a 


(2 a 1/3 y 


when \E — E c | < A 1/2 and \<p — <p c \ < 


^ 1/2 
cos E 


(ii) 


and zero otherwise, and A 1/2 = FWHM/2 is the WVR 
beam half- width. Note that as long as E c > A 1 / 2 ? ground 
pickup is avoided for this WVR beam pattern. By compar- 
ing the numerical results, the systematic error calculated 
using beam intensity distribution with sharp cutoffs (for 
the same FWHM) is about 20-30 percent s mal ler than 
for the Gaussian beams. (It was also found that neglect* 
ing beam spreading in the azimuthal direction underesti- 
mates the random error by less than 10 percent.) The 
radio telescope beam was approximated by an infinitely 
narrow pencil beam. Since the errors are proportional to 
tropospheric water vapor, all shown errors are for 1 g/cm 2 
of zenith water vapor column density. 


Before presenting numerical results, the WVR beam 
intensity distribution must be specified. The systematic 
error for an assumed Gaussian beam with 7.5-deg full 
width at half maximum (FWHM) has been calculated 
previously. 14 To simplify the computations, and since the 
aim of this article is to provide error estimates (rather than 
to tailor the errors to specific WVR beam designs), a beam 
is used whose cross section when viewed in the propaga- 
tion direction is a square with sharp cutoffs for the beam 
intensity. Specifically, for a beam centered at elevation E c 
and azi muth y> c , the WVR antenna temperature, T A (E c ) y 
is calculated by integrating the brightness temperature of 
the sky (Tb(E)) over the intensity distribution: 

T A (E e ) = [T B ] e ~ J J B e (E,<p) T b (E) dE cos E d<p 

( 10 ) 


A. Tropospheric Fluctuation-Induced Error for 
Brightness Centroid Pointing - 

The elevation E\, of the WVR beam brightness c entr oid 
is determined by requiring that the statistically averaged 
brightness temperature {Tb{E{)) at the centroid elevation 
Eh be equal to the WVR antenna temperature ( T A (E C )) 

(T B (E b )) = ( T A (E e )> (12) 

where E c is the elevation of the WVR beam geometrical 
center, and (...) designates the statistical average. Substi- 
tuting Tb(E) [Eq. (5)] into Eq. (10), neglecting the effect 
of ray bending and the Earth’s curvature, and using the 
optically thin troposphere approximation, the statistically 
averaged antenna temperature is 


where [. . ] c signifies the WVR beam average around 
(E Ci <p c ) } E and <p are elevation and azimuth angles, respec- 


12 J. Z. Wilcox, “The Error in the Estimated Path Delay due to WVR 
Antenna Beam Width: Beam Averaged Air Mass and Wet Tropo- 
sphere Fluctuations Effects,” JPL Interoffice Memorandum 335.6- 
92-004 (internal document), Jet Propulsion Laboratory, Pasadena, 
Califor nia, Janu ary 31, 1992. 

13 See Footnotes 4 and 5. 

14 S. Robinson, “A Simple Analytic Correction for WVR Beam 
Width,” JPL Interoffice Memorandum 335.4-530 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, July 23, 
1985. 


(T A (E e )) = T C + Tmc r, [A] e (13) 


where the WVR beam averaged air mass [A] c is given by 
the same integral expression as Eq. (10) except that Tb(E) 
is replaced by Ae = 1/sin E. By also using Eq. (5) for 
Ts(Eb) in Eq. (12), one obtains the result that in an op- 
tically thin troposphere, the brightness centroid coincides 
with the air mass centroid, 


1 

sin Eb 


/ 


BJ^p) 

sin E 


dE cos E dtp 


(14) 
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The calculated difference (the offset) between E c and 
E b is shown in Fig. 4. The offset is always positive (i.e., 
the brightness centroid is tilted from the beam’s geometri- 
cal center to a lower elevation), it increases with the beam 
width approximately as A 2 ^ 2 > an< ^ ^ has a ver y wide min- 
imum in the elevation range around E c ~ 52 deg. The 
minimum occurs because the systematic error, and hence 
the offset, depends on a nonvanishing second derivative of 
Ag versus E. Specifically, for E < 90 — Ax/ 2 » the offset 
E c - E b ~ (A 112/6) A” /A* e , where A f c and A” are the first 
and second derivatives of Ae — 1/ sin E versus E at E c . 
Since A" fA f c has a local minimum at 52 deg, the offset has 


also a local minimum at 52 deg. Note also that as E c ap- 
proaches 90 deg, the offset rapidly increases to Ax^/V^- 
For the current WVR Ax/ 2 < 4 deg, the offset is less than 
0.2 deg in the elevation range between approximately 20 
and 80 deg. Note that in an optically thick troposphere, 
the brightness centroid will differ somewhat from the air 
mass centroid. This is a consequence of the nonlinear de- 
pendence of Tb on the air mass, Eq. (5), in an optically 
thick atmosphere. 

The error in the estimated T B (E b ) is the square root of 
the variance 


<r$(E h ) = ((T A {E c )-T B (E b )) 2 ) = ( ( / T B (E) (B c (E,<p) - 6(E - E b ,<p - <p b )) dE cos E dip 


where S(E— E b ,<p—<p b ) is the Dirac delta function centered 
at (E b) (p b ). Equation (15) was evaluated by expressing 
T b (E) using Eq. (5), and then evaluating the correlations 
between the tropospheric opacities using the Kolmogorov 
turbulence model, as described in the paragraph following 
Eq. (7). The corresponding path delay error (£&)) 
was obtained by substituting <TT(E b ) into Eq. (2), where 
<r T (E b y& were identified with STbj's (j =1,2 designates 
20.6- and 31.4-GHz frequency channels, respectively) for 
the two WVR frequencies. Figure 5(a) plots the path de- 
lay error versus E c . Note that the error increases with 
decreasing elevations faster than L v . The error is plotted 
versus A 1/2 in Fig. 5(b). After a rapid rise near zero, the 
error increases sub linearly. At A 1/2 — 3 deg and E c = 30 
deg, 20 deg, and 10 deg, the errors are 0.06 cm, 0.1 cm, 
and 0.3 cm, respectively. Advanced WVR’s with narrow 
beamwidths were designed for observations at low eleva- 
tions. For Ax /2 — 1 deg, the errors are approximately 0.23 
cm and 0.5 cm at E c = 10 deg and 6 deg, respectively. 
These results qualitatively agree with errors calculated for 
copointed beams. 15 

The error shown in Fig. 4 refers to instantaneous mea- 
surements. However, astrometric data are averaged over 
time intervals on the order of 1 to 2 minutes, which tends 
to average out the fluctuations. Assuming that the radio 
telescope and the WVR observe simultaneously and con- 
tinuously during tint, the time averaged T B (E) is 

15 S. Keihm, “Finite Beam Effects on LOS Path Delay Decorrelar 
tion,” (internal document), Jet Propulsion Laboratory, Pasadena, 
California, March 22, 1990. 


T b (E) = -L / dt Tb(EJ) (16) 

tint JO 

Using T b (E) instead of T B (E) on the right-hand side 
of Eq. (15), the result is shown in Fig. 6. At ti nt < 
T 1 / 2 , the error decreases with a time constant T x / 2 — 
2 Ai/ 2 h v /v sin 2 !^, which is the time required for the 
moving troposphere to pass through the WVR beam cone 
(and thus erase the tropospheric differences between the 
beam averaged and line-of-sight opacities). The value of 
the time constant T x / 2 increases with decreasing elevation 
and increasing half-width. For h v = 2 km, v = 10 m/sec, 
Ax/a = 3 deg, and E c = 30 deg and 10 deg, T 1/2 is approx- 
imately 1.5 min and 12 min, respectively. For Ax/ 2 ~ 0.1 
deg, Tx /2 is 3 sec and 25 sec, respectively. Obviously, so 
that the time resolution is not degraded, WVR integration 
should never be longer than radio telescope integration. 

B. Tropospheric Fluctuation-Induced Error 
for Geometrical Center Pointing 

The systematic and random errors for a copointed radio 
telescope and a WVR antenna are calculated next. The 
usual argument for why WVR data should be associated 
with the direction of the beam’s the brightness centroid 
is that for a constant troposphere, (T^(E C )) = ( T B (E b )). 
If this were the only criterion, one could also correct the 
WVR data (i.e., Tj\(E c )) by subtracting from them the 
estimated value of the difference between (Ta(E c )} and 
(T b (E c )) and identify this “corrected” data as the actual 
value of the brightness temperature T B (E C ) in the direc- 
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tion of the beam’s geometrical center. For a copointed 
radio telescope and a WVR antenna, the (systematic) dif- 
ference AT A (E c ) y is determined as 




A ?/2 
6 sin E c 



_ 2 -) 
tan 2 E c J 


( 18 ) 


AT a (E c ) = (T a (E c ) - T b (E c )} = T mc t z A A e (17) 


where r z is the total (wet and dry) zenith opacity, and 
AA C = [A\c — A{E C ) is the difference between the beam 
averaged and geometrical center air mass. Note that 
{T A (E c )) is bigger than ( T B (E C )), 

The value of A T A (E C ) at 20.6 GHz (and the corre- 
sponding path delay error) is shown in Fig. 7. The error 
increases with beamwidth approximately as A 2 / 2 . Note 
that for the radiation pattern with a sharp cutoff for in- 
tensity distribution, this quadratic dependence on A 1/2 
can be derived analytically 


where the approximate equality on the right-hand side has 
been obtained for narrow line widths, A 1/2 << E c . The 
more important feature to notice in Fig. 7 is that the e rror 
increases with decreasing E c more rapidly than L V1 namely 
that A L v oc AT a (E c ) oc AA c oc (1 + 2/ tan 2 .E c )7sin E t . 
This is the same type of increase as for the random error 
calculated in the preceeding paragraphs, i.e., the errors 
caused by WVR beam averaging increase with decreas- 
ing elevation more rapidly than a simple scale error. At 
Ax /2 — 3 deg and E c = 30 deg, 20 deg, and 10 deg, the 
systematic path delay error is approximately 0.36 mm, 1.2 
mm, and 10 mm, respectively. The error decreases with 
decreasing beam width. At Ai/ 2 1 deg and E c = 10 
deg and 6 deg, the error is 2 mm and 6 mm, respectively. 
(At Ax /2 — 04 deg, the error is 0.1 mm and 0.06 mm, 
respectively.) 


AA C = — 
c 2 A 


1/2 


tan(£ c + A 1/2 )/2 
tan (E c — A 1/2 )/2 


Correcting the measured T A (E e ) by A T A (E C ), the sto- 
chastic error in the inferred Tb{E c ) is determined from 


I 


4(e c ) = ((t a (e c )-at a (e c )-t b (e c )) 2 } 


C Tb(E ) - ( T b (E ))) (B C (E, <p) - 6{E -E c ,<p- <p c )) dE cos E d<p 


(19) 


where 6(E — E Ci <p — <p e ) is the Dirac delta function cen- 
tered at (E C} <p c ), Equation (19) was evaluated by using 
the same procedure as Eq. (15). The corresponding path 
delay error (ai )V (E c )) is shown in Fig. 8(a) versus E c and 
in Figure 8(b) versus Ax/ 2 - Similarly as for the system- 
atic error [and for the stochastic error for £HeT)rightness 
centroid, (£&)], <tl ^(Ec) increases with decreasing E c 
more rapidly than L v . Note that when the systematic error 
is smaller than (Ti tV (Ef > ) (such as for Ax/ 2 j< 1 deg at a 6- 
deg elevation, or for Ax / 2 < 1.5 deg at a 10- deg elevation), 
<rL,v(Ec) is about the same as <ri iV (Eb ); whereas when the 
systematic error is bigger than <7£ >t ,(£&), ai tV (E c ) looks 
more like the systematic error (which increases oc A 2 / 2 and 
can become very large). That is, while the copointing will 
significantly increase the stochastic error for wide WVR 
beams, the increase will be small for narrow beams. It 
has also been found that when <tl }V (E c ) ~ <ri )V (Eb) (the 
narrow beam case), the (Tl,v{E c ) for the integrated sig- 
nal [as in Eq. (16)] decreases with tint at about the same 
rate as does (?L } v{Eb)\ whereas when <tl )V (E c ) > (Ti iV (E{) 


(the wide beam case), the decrease is significantly slower. 16 
Thus, when the systematic error is less than <tl,v{Ei>) (the 
narrow beam case, Ax/ 2 < 1 deg for all E > 6 deg), the 
copointing will introduce a negligible error into the esti- 
mated path delays using WVR’s. However, for beam sizes 
greater than 1 deg, the WVR’s should be pointed at a 
slightly higher elevation than the radio telescope. 

V. Discussion and Recommendations 

The main goal of this article is to investigate how tro- 
pospheric dynamics affect the ability of realistic WVR’s to 
track tropospheric fluctuations. Two effects were studied 
in detail: errors in WVR instrument gain cal ibration from 
tip curves and errors in the estimated brightness temper- 
ature caused by WVR beamwidth averaging. The errors 
can be used to derive WVR gain stability requirements 


16 See Footnote 12. 
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and WVR antenna beamwidth that would make it possi- 
ble to reduce the path delay error to the 1-mm level in the 
elevation range from zenith to 6 deg. 

The minimum error induced by tropospheric fluctua- 
tions in a single gain estimate (per 1 g/cm 2 of zenith 
water vapor column density during WVR calibration) is 
approximately 0.26 percent. That error causes two types 
of errors in the estimated path delay. The first error, ap- 
proximately 3 mm, is independent of path delay. Provided 
that the WVR gain remains constant, this is a bias error 
that can be removed by differencing between VLBI obser- 
vations (biases have no effect on delay rates used for gravi- 
tational wave searches). The second error, 6L V ~ L v (Tg/g , 
is a scale error. Figure 9 shows the scale error for a sin- 
gle gain estimate 6L V (mm) ~ 0.16 N ViCa i^v,ob»/ s ' m & 
(where N VtCai and N Vf0 bs come from water vapor con- 
tent during the WVR calibration and radiometric obser- 
vation) as a function of N VtCa iN Vt0 bs and elevation ( E ) 
of the radiometric observation. For example, assuming 
that N VtC ai = 1 g/cm 2 (corresponding to L VfZ = 6 cm) 
and N v>obt = 2 g/cm 2 , SL V exceeds 3 mm (i.e., it ex- 
ceeds the Lv-independent error) when E < 6 deg. When 
N V} cai = N Vf obs ~ 2 g/cm 2 , SL V exceeds 6.2 mm at E = 6 
deg. To achieve the desired 1-mm path delay accuracy, ei- 
ther the effect of SL V on astrometric estimates or Sg itself 
must be reduced. The success of any approach to obtain- 
ing accurate astrometric estimates depends on the stability 
of the WVR gain. 

For a stable WVR gain, the scale error is systematic. 
Preliminary results of attempts to reduce the effect of sys- 
tematic scale errors on astrometric estimates with VLBI 
data analysis appear to be promising, although more work 
is needed to ascertain quantitative results. 17 The gain 
error can be reduced by using an alternate gain calibra- 
tion technique (such as two absolute reference load cal- 
ibrations), or, assuming that the error induced by tro- 
pospheric fluctuations is the dominant error source, by 
tip curve repetition. Uncertainties in alternate calibra- 
tion methods have so far prevented circumvention of tip 
curves. To reduce the error by tip curve repetition, the 
WVR gain must be sufficiently stable. For example, to 
achieve a 1-mm path delay accuracy at a 6-deg elevation 
when N Vt obs = 2 g/cm 2 , the gain error must be less than 
0.08 percent. To reduce the gain error to 0.08 percent 
when N VtCa i = 2 g/cm 2 , the tip curve must be repeated at 
least (0.52/0.08) 2 = 40 times. Section III showed that suc- 
cessive gain estimates become decor related within a typ- 


17 R. Linfield, Tracking Systems and Applications Section, personal 
communication, Jet Propulsion Laboratory, Pasadena, California, 
1991. 


ical Tcorr of approximately 7 min. Therefore, if one at- 
tempts to reduce the gain error by tip curve repetition, the 
tip curves should be separated by a time interval greater 
than 7 min, and the WVR gain should change by no more 
than 0.08 percent over at least 40 T corr ss 5 hr. 18 (Be- 
cause of the effect of other error sources and the possi- 
bility of humidity higher than 2 g/cm 2 , the recommended 
WVR stability is 0.04 percent over the 5-hour period.) In 
the time interval between the WVR calibration and the 
radio metric observations (and during the observations) 
the gain will still have to be updated, e.g., by comparing 
the number of counts for the reference load. 19 ’ 20 

The error in the estimated brightness temperature due 
to WVR beam averaging of tropospheric fluctuations was 
found to be smaller when the direction of the radio source 
coincided with the WVR beam brightness centroid than 
with the beam geometrical center. The error increases 
with tropospheric water vapor content and beamwidth. 
More important, however, is that the error increases with 
decreasing elevation faster than L v . Since low-elevation 
data will be weighted more heavily than they should be, 
this will affect astrometric estimates. Advanced WVR de- 
signs have been suggested to reduce the error and avoid 
ground pickup by implementing narrow beams. Figure 10 
shows the errors in the two-dimensional space of beam 
half-widths (Ai/ 2 ) and elevations for 1 g/cm 2 of zenith 
water vapor column density. For each curve, the error is 
less than the cutoff error for all Ai/^s and & s below and 
to the right of the curve. For example, for the error to 
be less than 1 mm at all E > 10 deg, A x / 2 should be 
< 0.1 deg. For more humid weather, the errors will be 
higher (and the beamwidth requirement more stringent), 
proportional to zenith water vapor. Because of various ap- 
proximations involved in deriving Fig. 10 (the WVR beam 
radiation pattern with sharp cutoffs for intensity distribu- 
tion, an infinitely narrow pencil beam for the radio tele- 
scope, and an optically thin troposphere), the guidelines 
are approximate (the guidelines can easily be quantified by 
applying the methods described in this article to specific 
beam shapes). 

Signal integration reduces the fluctuation-induced er- 
ror from its instantaneous value with a time constant 
Tij 2 — 2 Ax /2 hv/ v sin 2 E. Because of its dependence 
on beam width and elevation, the effect of signal integra- 


ls The measured gain of the current J and D series WVR s drift 
at a rate that causes an approximate 0,3-percent gain change in 
1000 sec (Footnotes 4 and 5). Therefore, the stability of these 
WVR’s should be improved by a factor of at least 60, from g/g ~ 
3 X 10” 6 sec” 1 to i/^5X 10 -8 sec -1 . 

19 See Footnote 2. 

20 See Footnote 4. 
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tion has been neglected in Fig. 10. For example, for a 
WVR beamwidth A X j 2 = 0.1 deg (FWHM = 0.2 deg) and 
E = 6 deg, 7|/2 is 1 min, and the path delay error when 
N Vi obs = 2 g/cm 2 is about 4 mm. Hence, integrating the 
WVR signal over a 2-min period (which is a typical VLBI 
integration time) will reduce the fluctuation error to about 
1 mm (which is the desired accuracy for the estimated path 
delays). It has heen suggested that copointing a radio tele- 
scope and a WVR antenna in the same direction would 
simplify WVR antenna steering. The copointing will in- 
troduce a systematic error and increase the random error 
in the estimated path delay. For E > 6 deg, these addi- 
tional errors will be smaller than the random error for the 
brightness centroid pointing for all beam sizes A X j 2 < 1 
deg. For A X f 2 > 1 deg, the additional errors will increase 
oc A 2 /2 , and, in addition, the required WVR signal inte- 
gration time to average out the fluctuation-induced error 
becomes longer 21 than the T 1/2 for the brightness centroid 
pointing (and longer than the VLBI integration time of 
about 2 min). Therefore, the ability for the wide beam 
WVR’s to be pointed at a slightly higher elevation than 
the radio telescope is important. 


21 See Footnote 12. 


From the brief discussion of various error sources in Sec- 
tion II, the biggest error in the estimated path delay is at 
the present time due to inaccurate modeling of the absorp- 
tivity of water vapor and uncertainties in the distribution 
of atmospheric parameters along the observed lines of sight 
(atmospheric noise). These uncertainties cause scale errors 
(see Fig. 1) in the estimated path delay: a systematic er- 
ror of about 10 percent due to the error in the absorptivity 
model and between 2 and 4 percent random error due to 
atmospheric noise. To satisfy the 1-mm path delay ac- 
curacy requirement at a 6-deg elevation, the atmospheric 
noise must be reduced to the 0. 17-percent level by, for ex- 
ample, custom tailoring the retrieval algorithm coefficients 
to specific sites and a set of observing conditions. 

The accuracy of the present absorptivity models could 
be improved by better modeling and model calibration, 
using for example, a comparison of WVR and radiosonde 
data, direct measurements or estimates using interferomet- 
ric data reduction of atmospheric path delays, or mea- 
surement of water vapor absorptivity in a laboratory- 
controlled environment. Investigation of some of these pos- 
sibilities, including that of developing mathematical meth- 
ods to filter out the effect of the systematic scale error 
during VLBI data reduction, is part of an ongoing effort. 
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Fig. 1. The effect of stochastic (e.g., due to thermai noise and 
WVR gain variation) and scale (e.g., from absorption coefficient 
and WVR calibration) errors on path-delay estimates. The scale 
errors plotted are per 1 g/cm 2 of water vapor zenith column den- 
sity (or, equivalently, per 6 cm of zenith path delay). The goal is 
1-mm path-delay calibration accuracy in most weather conditions. 
(This figure assumes no reduction in the effects of a scale error 
from parameter estimation.) 



Fig. 2. The standard deviation of a single gain estimate from 
tip curves versus minimum tip elevation. The heavy and broken 
line curves were calculated using the actual and unit observable 
covariance-variance matrix W” 1 in the least-squares analysis. 
The path delay shown Is the bias error A L v . 



Fig. 3. Temporal development of the square root of the covariance 
of successive gain estimates versus time, 7, elapsed between 
tip curves. The single gain estimates decorreiate within the time 
Tcorr = hy/(v tan E m f n ). The broken line curve connects co- 
variance values at 7 = Tcorr/ 2. 
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PATH DELAY BIAS ERROR, mm 









Fig. 4. The offset (A E c ) between elevations of the WVR beam ge- 
ometrical center and brightness centroid. For small beam widths, 
A E c oc A 2 j 2 * Note that as E c approaches 90 deg, the offset 
Increases to A 1 / 2 / V3- 




Fig. 6. Path delay stochastic error for the brightness centroid 
pointing versus the signal Integration time. For a symmetrical 
radiation pattern, the error depends only very little on the wind 
direction (for different wind directions, the calculated errors dif- 
fered by less than 15 percent). 


Fig. 5. Path delay stochastic error for the brightness centroid 
pointing due to WVR beam averaging of tropospheric fluctuations 
(Instantaneous error, per 1 g/cm 2 of zenith water vapor column 
density): (a) versus elevation and (b) versus beam half-width. 
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Fig. 7. The systematic difference A T& between the antenna tem- 
perature and the brightness temperature at 20.6 GHz (and the cor- 
responding path delay error) In the beam geometrical center due 
to WVR beam width averaging of air mass (per 1 g/cm 2 of zenith 
water vapor column density): (a) versus elevation and (b) versus 
beam half-width. 



BEAM HALFWIDTH, deg 


Fig. 8. Path delay stochastic error for the geometrical center 
pointing due to WVR beam averaging of tropospheric fluctuations 
(per 1 g/cm 2 of zenith water vapor column density): (a) versus el- 
evation and (b) versus beam half-width. ' 


48 







Fig. 9. Curves of path delay scale errors (5 L v ) for a single gain 
estimate. For each curve, the error Is less than the cutoff error 
above and to the left of the curve. The values of Ny^and N v obs are 
zenith water vapor column densities during WVft calibration and 
radiometric observation, respectively. 



Fig. 10. Curves of tropospheric fluctuation-induced path delay 
error [g l v (E^) for the brightness centroid (per 1 g/cm 2 of zenith 
water vapor column density) due to WVR beam averaging. In more 
humid weather, the errors scale with humidity. 
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Appendix A 

The Method of Least Squares 


For brevity of notation, the voltage recorded by a WVR 
pointed in a direction E 9 is designated K*, and the line-of- 
sight opacity is designated r,*. The duration of a single tip 
sequence is less than the time it takes to decorrelate the 
tropospheric inhomogeneities. The ith tip curve voltage, 
Vi , is then modeled as 

Vi = g{T c e- T < -T M {l-e- T ')-T rtS ) 

— 9 T i T\fc — g T rc (A- 1) 

where Trc = T rt j - Tc, Tmc = Tm - T Ci and where 
the linearized Eq. (A-l) is a good approximation to the 
full radiation transport equation for most (r,- < 0.5) tro- 
pospheres of interest. 

In tip curve analyses, the WVR recorded data are fit 
by assuming a temporally constant and spatially homoge- 
neous stratified troposphere. For i = 1,... N data, where 
N is the number of tip curve elevations, the V^s represent 
a set of N equations for two solve- for parameters, g and f 2 . 
Mapping the r f -’s to zenith by using r* = f, , A* , where r z 
is the estimate for the zenith opacity r Z) and the air mass 
A{ = 1 / sin Ej } the equations are solved by the method of 
least squares [6], as follows. 

Defining solve-for parameters and observable column 
vectors X = [g r z , g] and V = [Vi,... Vjv], respectively, 
the design matrix A (A has dimensions N x 2) in V = A X 
is 

Tmc A\ — Trc \ 

TmcA2 —Trc 

: : ( A "2) 

Ta/cAw -Tr C ) 

Assuming that the errors in V*’s have zero means and 
a variance-covariance matrix W~* 9 minimization of the 


quadratic form ((V - AX) T W (V - AX)) yields the fol- j 

lowing estimates for g and r z : i 

X = (A t WA)- 1 A t W V (A-3) 1 

where X is the column vector X — [gf Zy y], and the super- I 

script T designates the transpose matrix. By substituting I 

Eq. (A-l) into Eq. (A-3), one can easily verify that the 1 

statistically averaged estimated gain is equal to the WVR i 

gain (i.e., {g) = g) } as it should be. The estimated gain 
standard deviation is the square root of a ? = cov(^,^), 
which is the matrix element (o^) 2,2 of 

4 S Exp {(X - X)(x - Xf} 

= B~ 1 A t W cov(F,K t ) W A B~ l (A-4) = 

where Exp designates the expectation value, cov(V, V T ) 
is the actual observable covariance-variance matrix, and 
B = A t WA. Equations (A-3) and (A-4) yield g , r Zf and - 
c 7x for given Wij and cov(V, V T ). 

In practice, X and <j\ can be derived either by us- 
ing some assumed W™ 1 (the so-called consider analysis) 

[1], or by setting W~ l equal to the observable variance- 

covariance matrix cov(V, V T ). The most common (and 

simplest) form of the assumed W" 1 is the unit matrix = 

(W~ 1 )ij = 6 if j , where Sij is the Kronecker delta. Taking 

the unit W~ l corresponds to assuming that the observable 

errors are uncorrelated with equal variances and reduces 

the minimization procedure to the minimization of the sum 

of squares. The latter, i.e., setting (W~ 

minimizes the variance-covariance matrix cr\ as — 

(A T WA)-K 

The results of these calculations are discussed in the 
main text. Evaluation ofcov(V*, V,) using the Kolmogorov 
turbulence model is described in Appendix B. 
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Appendix B 

Evaluation of the Correlation Functions Using the 
Kolmogorov Turbulence Model 


To evaluate correlations between the simulated data, 
the voltage Vi associated with the sky brightness in 
the elevation direction E{ is first related to the line-of-sight 
opacity r< by using Eq. (A-l). This expresses correlations 
between VJ’s in terms of correlations between tv’s. Next, 
neglecting fluctuations in the dry component of r», cor- 
relations between wet opacities are expressed in terms of 
correlations of wet refractivity by expressing the wet opac- 
ity (r V) i) as the line-of-sight integral 


,. = / 


h v j sin Ei 


l>h v lainE t 

0.(3) dn ~ / \{?i)dri 

L V> Z Jo 

(B-l) 


order of integration and ensemble averaging, and then set- 
ting dr { = Ai dz and d r, = Aj dz ' , Eq. (B-2) yields 



where the variance <r 2 of the wet refractivity fluctuation 
is independent of spatial coordinates and is obtained by 
letting the distance R go to infinity in 


where h v is the height of the tropospheric slab; ot v is the 
water vapor absorptivity per unit length; r VfZ and L VfZ 
are the average values of wet opacity and path delay at 
zenith, respectively; and x(^») 19 index of refraction — 
1. The correlations between x(^)’ s are evaluated by us- 
ing the Kolmogorov turbulence structure function, Eq. (8) 
of the main text. In what follows, the above described 
procedure of evaluating the correlations is exemplified in 
the evaluation of the observable variance-covariance ma- 
trix W~ l . 

By using Eq. (A-l), the matrix element W~j is written 
as 


W~> = cav(Vi,Vj) = g 2 Tif C cov(T v>i ,T vJ ) 

= 9*1% c T *j) ~ ( B_ 2) 


= (<x’> - «’) 


= CO ) _ N% C 2 lV 3 

2 “ 2 

(B-5) 


where the last expression on the right-hand side has been 
obtained by evaluating the asymptotic D X (R = oo) by us- 
ing Eq. (8) of the main text, and where N v is the water 
vapor column density at zenith in g/cm 2 , and L 9 is the tro- 
pospheric turbulence saturation length. The reason why 
D x ( oo) should converge as R becomes very large has been 
discussed in [1]. 

The covariance of successive gain estimates and the 
effect of signal integration on beam averaging of tropo- 
spheric fluctuations involve evaluation of correlations of 
type (Ti(t)Tj(t + T)). By using the “frozen” troposphere 
model [1], these correlations were evaluated using the ex- 
pression 


where (...) signifies statistical ensemble averaging. 

By using Eq. (B-l) and the expression (Eq. A.3 of [1]) 

(XiXj) = <x 2 > - \ D x {n - 3) (B-3) 

where x« = x(3) and D x {?i - fj) = ((x« - Xj) 2 ) * s l ^e 
refractivity spatial structure function by interchanging the 


(x(ri,t)x(rj,t + T)) = (x 2 ) - \ D xi^ ~ +^ T ) ( B ‘ 6 ) 

where the structure function 

D x (fi -fj+v T) = ((x(3,0-x(3» < + 7 ’)) 2 ) 

is the same as Eq. (8) of the main text. 
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